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ABSTRACT 

A  two-dimensional  transient  computer  program  developed  for  fluidization 
was  used  to  study  the  uniformity  of  powders  injected  into  rectangular 
cavities.  Most  configurations  gave  very  non-uniform  powder  distributions. 
Simulation  of  a  leak  present  in  the  experimental  bag  experiments  produced 
a  jet  expansion  not  unlike  that  observed  in  photographic  records.  Simula¬ 
tions  with  5,  30  and  200  micron  powders  showed  that  the  jet  expansion  was 
most  pronounced  with  5-micron  powders,  and  that  the  powder  distribution 
was  the  most  uniform  of  the  configurations  studied. 

A  comparison  of  experimental  concentration  contours  obtained  from 
the  high-speed  photographic  records  of  powder  injection  showed  a  reason¬ 
ably  good  agreement  with  the  theory.  The  agreement  is  considered  good, 
because  the  theory  has  no  fitted  parameters.  It  is  predictive.  Discrep¬ 
ancies  between  the  theoretical  and  experimental  results  may  be  qualita¬ 
tively  accounted  for  by  bag  expansion  and  a  reduced  pressure  in  the  bag 


due  to  leaks. 
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INTRODUCTION 

IITRI  bag  experiments  which  involve  injection  of  powders  into 
partially-collapsed  plastic  bags  were  simulated  using  a  transient  two- 
dimensional  two-phase  computer  program  (Gidaspow  and  Ettehadieh,  1983). 
The  objective  was  to  determine  the  uniformity  of  powders  for  various 
injection  modes.  The  bags  were  approximated  by  a  two-dimensional  slice 
with  rigid  walls.  Since  the  rigid  walls  may  not  properly  simulate  a  bag 
that  inflates  upon  the  injection  of  the  powder-air  streams,  leaky  walls 
were  also  used  in  modeling,  to  give  a  constant  pressure  in  the  two  dimen¬ 
sional  cavity. 

MATHEMATICAL  MODEL 

The  particle  injection  model  uses  the  principles  of  conservation  of 
mass  and  momentum  for  each  phase.  It  is  a  two-fluid  model  of  the  type 
reviewed  by  Lyczkowski,  Gidaspow  and  Solbrig  (1982)  and  by  Gidaspow, 
Rasouli  and  Shin  (1983) .  Although  there  are  no  studies  known  to  the 
writers  dealing  with  hydrodynamic  modeling  of  powder  injections  into 
containers,  the  modeling  dealing  with  propellant  motion  by  Koo  and  Kuo 
(1977)  and  by  Gough  (1979)  is  similar  to  the  present  approach.  The 
equations  used  for  powder  Injection  modeling  are  those  developed  by 
Gidaspow  and  Ettehadieh  (1983)  for  modeling  fluidization. 

Table  I  shows  the  continuity  and  the  separate  phase  momentum 

equations  for  two-dimensional  transient  two-phase  flow.  There  are  six 

nonlinear  coupled  partial  differential  equations  for  six  dependent 

variables.  The  variables  to  be  computed  are  the  void  fraction,  e,  the 

pressure,  P,  the  gas  velocity  components  U  ,  V  ,and  the  solids  velocity 
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components  U  and  V  ,  in  the  x  and  y  directions,  respectively.  The 
equations  are  written  in  a  form  found  in  the  K-FIX  computer  code 
(Rivard  and  Torrey,  1977)  for  gas-liquid  flow. 


Table  I  Hydrodynamic  Model 
CONTINUITY  EQUATIONS 


Gas  Phase 


at  ^pge^  +  3x  pgUg^  +  9y  pgV  ~  0 


Solid  Phase 


at 

MOMENTUM  EQUATIONS 


hr  [Pjl-e)l  +  y  [PSUS<1-01  +  h  [Psd-e)Vs]  -  0 


3y  lKV 


Gas  Momentum  in  x  Direction 
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Solids  Momentum  in  x  Direction 
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Gas  Momentum  in  y  Direction 
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Fluid  Particle  Drag 

The  major  empirical  input  in  this  model  is  the  fluid-particle  friction  co 
efficient,  6*  It  is  obtained  from  standard  correlations  as  explained  in 
detail  by  Gidaspov  and  Ettehadieh  (1983) . 

For  porosities  less  than  0.8  the  pressure  drop  due  to  friction  be¬ 
tween  gas  and  solids  can  be  described  by  the  Ergun  equation.  The  friction 


coefficient  in  this  porosity  range  becomes 
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For  porosities  greater  than  0.8  Stokes'  law  corrected  for  high  flow 
and  packing  leads  to  the  following  expressions  for  the  friction  coefficient : 
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where  C^,  the  drag  coefficient  in  the  y  direction,  is  related  to  the 
Reynolds  number  by  means  of  the  well  known  (Soo,  1967)  relations  given 
below. 
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In  equation  (2)  f(e)  shows  the  effect  due  to  the  presence  of  other  particles 
in  the  fluid  and  acts  as  a  correction  to  the  usual  Stokes  law  for  free 


fall  of  a  single  particle.  We  use 


In  the  present  study  this  correction  is  near  one  and  therefore  contributes 
very  little  to  the  results.  The  expression  for  the  friction  coefficient 
in  the  x  direction  is  the  same  as  that  in  the  axial  direction,  y. 

Solids  Stress 

In  a  general  formulation  the  solids  momentum  equations  would  contain 
solids  stress  terms  that  are  a  function  of  porosity,  pressure  and  the 
displacement  tensors  of  solids  velocity,  gas  velocity  and  relative  velocity 
Such  a  general  formulation  with  proper  values  of  material  constants  does 
not  exist  today.  However,  it  is  necessary  to  use  the  normal  component  of 
the  solids  stress  to  prevent  the  particles  from  reaching  impossibly  low 
values  of  void  fraction.  Measurements  of  such  a  stress  for  settling 
showed  it  to  be  quite  small  compared  with  the  hydrostatic  pressure.  In 
view  of  the  previous  theory  and  measurements,  the  constitutive  equation 
for  the  normal  component  of  the  stress,  T  is 

x  «  T(e) 


Then  by  chain  rule, 
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Let  the  modulus  of  elasticity  or  particle  to  particle  interaction  co¬ 
efficient  be 

G(e)  -  ||  (8) 


In  this  study  this  modulus  is  fitted  to  the  experimental  data  for  fluidiza¬ 
tion  (Gidaspow  and  Ettehadieh,  1983)  by  the  equation 


-G(e)  -  io("8,76e  +  5*43>  ,  N/m2  (9) 
This  term  becomes  of  numerical  significance  only  when  the  void  fractions 
go  below  the  minimum  fluidization  void  fraction.  Although  this  does  not 
happen  in  this  study  this  term  helps  to  make  the  system  numerically 
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stable,  because  It  converts  the  Imaginary  characteristics  into  real  values. 

Fanucci,  Ness  and  Yen  (1979)  have  obtained  the  characteristic  direc¬ 
tions  for  one  dimensional  incompressible  flow  for  the  set  of  equations 
shown  in  Table  I.  Their  paper  shows  one  of  their  characteristics,  C  to 
be 


The  above  relation  shows  that  as  the  void  fraction  goes  to  zero  -G(e) 
becomes  large  and  makes  the  characteristics  real.  In  practice  we  find 
that  this  extra  stress  permits  us  to  continue  calculations  in  time  with¬ 
out  a  loss  of  material  balance.  Without  this  term  the  K-FIX  code  (Rivard 
and  Torrey,  1977)  gave  results  that  showed  a  loss  of  mass  after  about  0.7 
seconds  of  operation  even  with  a  very  strict  convergence  criterion. 


Boundary  and  Initial  Conditions 

Initially  there  was  no  powder  in  the  bag.  Thus  the  powder  volume 
fraction  is  set  to  zero.  The  powder  of  a  prescribed  volume  fraction 
enters  the  system  at  prescribed  solid  velocities.  Since  these  velocities 
are  high  and  the  inlet  gas  flow  rates  are  not  known,  it  is  assumed  that 
there  exists  no  slip  between  the  powder  and  the  air  entering  the  bag. 

Thu6  the  gas  velocity  entering  the  system  is  assumed  to  be  equal  to  the 
velocity  of  the  fine  powder.  This  assumption  may  not  be  exactly  correct 
«nH  may  lead  to  the  injection  of  more  air  than  is  true  in  the  IITRI  bag 
experiments.  A  knowledge  of  the  pressure  in  the  bag  before  explosion 
would  have  been  helpful  in  checking  this  assumption. 

The  boundary  conditions  (B.C.)  are  as  follows  (See  Figure  1): 


B.C.l  At  x-0,  the  prescribed  variables  are: 

At  the  inlet  port,  the  velocities  were  constant  or  as  shown  below. 


U  «  U  «  40  +  10.5  sin(20nt)  t<t 
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where  tQ  Is  the  time  at  which  the  powder  flow  was  terminated 

V  -  V  -  0 
8  s 

e  -  0.998 


Elsewhere  at  the  walls, 
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B.C.2 


For  configurations  such  as  "d"  (symmetrical) • 


At  x  *  200  cm 
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For  configurations  such  as  "f " , 
At  x  *  400  cm 
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B.C.3 


At  y  -  0 
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B.C.4 


At  y  ■  88  cm 


For  configurations  without  any  leak. 
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For  configurations  with  a  leak 

P  ■  atmospheric 

V  -  0 

s 

Numerics  and  Computational  Time 

The  two-phase  flow  equations  describing  the  dynamics  of  the  powder¬ 
laden  gas  were  solved  using  a  modified  K-FIX  program  (Rivard  and  Torrey, 
1977).  The  K-FIX  program  solves  these  equations  using  a  semi-implicit 
technique.  The  continuity  equations  are  differenced  fully  implicity, 
whereas  the  pressure-term  and  the  drag-term  are  the  only  implicit  terms  in 
the  differenced  momentum  equations.  At  any  time  step  these  difference 
equations  are  solved  iteratively,  to  obtain  the  updated  values  of  the 
dependent  variables. 

The  computations  were  carried  out,  typically,  over  a  two  dimension¬ 
al  region  consisting  of  150  computational  cells.  The  horizontal  direc¬ 
tion  was  discretized  into  ten  parts  of  40  cm  each  and  the  vertical 
direction  was  discretized  into  15  parts  of  a  height  of  6.1  cm.  The 
simulation  of  the  powder  flow  inside  the  bag  was  typically  carried  out 
for  two  seconds  of  real  time.  On  a  commercial  CRAY-I  computer,  about 
two  minutes  of  CPU  time  was  expended  for  each  of  those  computations. 

The  density  plots  shown  in  Figures  (2)  and  (8)  are  generated  from 
the  powder  volume  fraction  data.  The  number  of  dots  in  any  computational 
cell  is  proportional  to  the  square  of  the  powder  volume  fraction  in  that 
cell.  In  a  visual  observation  of  the  actual  process  such  a  nonlinear 
functional  relationship  between  the  powder  volume  fraction  and  the  total 
projected  area  is  believed  to  exist. 


COMPUTER  RESULTS 


Figure  1  depicts  the  various  bag  configurations  for  which  computer 
runs  were  made  together  with  the  inlet  conditions  that  were  used. 

Initially  the  bags  contained  only  air  at  one  atmosphere.  Figure  2  shows 
the  powder  concentration  at  various  times  for  configuration  "a".  The 
powder  moved  into  the  bag  with  a  slug-type  motion,  was  reflected  from  the 
wall  opposite  the  inlet  and  then  was  returned  to  the  inlet  wall.  Figure  2 
shows  that  there  exists  a  central  void  in  the  cavity  containing  no  powder. 
This  void  persisted  even  after  continued  injection  of  powder  beyond  the 
0.8  sec  shown  in  the  figure. 

Cases  b  and  d'  with  opposing  jets  were  run  by  assuming  a  central 
symmetry.  The  computer  program  was  not  programmed  to  handle  the  non- 
symmetrical  case  C,  because  the  computations  involve  a  sweep  from  one 
direction.  This  case  would  require  reprogramming. 

EFFECT  OF  LEAK,  SINUSOIDAL  INPUT  AND  MIXING 

The  introductory  runs  had  shown  the  powder  motion  to  be  of  an  almost 
slug  type.  It  was  suspected  that  this  was  in  part  due  to  the  build  up 
of  pressure  in  the  cavity  which  prevented  a  large  lateral  particle  motion. 
For  example,  for  configuration  d  the  pressure  rose  to  4.8  atmospheres 
after  0.7  seconds  of  powder  injection.  In  the  IITRI  bag  experiment  the 
pressure  is  clearly  not  that  high  because  of  bag  expansion  and  leaks. 

To  simulate  this  situation  properly  a  leak  along  the  upper  boundary  was 
allowed  in  configuration  "e".  The  leak  was  such  that  the  pressure  in  the 
bag  stayed  at  essentially  one  atmosphere. 

To  further  correct  for  greater  jet  expansion  in  the  lateral  direc¬ 
tion,  a  sinusoidal  variation  of  the  inlet  velocity  was  taken.  Figures  3a 
and  3b  show  that  the  sinusoidal  variation  had  a  negligible  effect  on  the 


jet  expansion.  However,  a  comparison  of  Figures  3a  and  3b  to  Figures  3c 
and  3d  shows  that  the  leak  had  a  very  large  effect  on  the  lateral  jet 
expansion.  Nevertheless,  the  run  with  the  leak  shows  that  the  powder 
continued  to  move  horizontally  as  well  as  vertically  in  the  direction  of 
the  leak;  a  more  complex  behavior  was  obtained  due  to  the  lack  of  res¬ 
training  pressure  that  had  been  present  in  the  runs  with  no  leaks. 

Figure  4  shows  the  gas  velocities  for  configuration  d,  with  no  leak 
The  gas  flow  reaches  an  almost  steady  circulatory  motion  after  about  0.1 
seconds  of  powder  injection.  The  solid,  on  the  other  hand,  lags  c r 
siderably  behind  the  gas  motion;  as  indicated  previously  in  Figure  a 
0.2  seconds  there  was  no  solid  present  along  the  top  and  bottom  w;  s. 

Figure  5  shows  the  gas  velocities  with  a  leak  for  configuration  e.  .»e 

gas  velocities  on  the  top  face  are  considerably  modified  by  the  leak  in 
the  expected  manner,  that  is,  flowing  out  of  the  cavity.  The  solid 
partially  moved  up,  but  did  not  follow  the  gas  stream  lines. 

Figure  6  shows  the  concentration  profiles  with  and  without  leak  at 
0.1  seconds  at  various  positions  in  the  cavity.  Generally,  the  solid 
moved  toward  the  leak,  but  there  exists  a  wave  motion.  For  example,  in 
Figure  6b  the  solid  moved  away  from  the  leak  and  tended  to  accumulate 
near  the  boundary  opposite  the  leak.  Figure  7  shows  the  concentration 
profiles  after  0.2  seconds  of  powder  injection.  The  accumulation  of 
solids  is  more  pronounced  due  to  the  wave  motion. 

In  runs  "d"  to  "f"  the  powder  was  terminated  at  0.82  seconds  and 
at  0.53  seconds,  but  the  powder  was  allowed  to  mix  due  to  inertia  of 
the  particles  rotating  in  the  cavity.  An  examination  of  the  computer 
output  of  solids  volume  fractions  showed  that  a  great  uniformity  of 
particles  was  produced  by  this  mode  of  operation.  However,  there  was 


a  considerable  settling  and  accumulation  of  powder  at  the  bottom  of  the 
cavity.  See  Table  II. 

PARTICLE  SIZE  EFFECT 

The  effect  of  particle  size  was  studied  in  configurations  "g"  to 
"i".  Figure  8  shows  the  particle  concentrations  at  0.1  and  0.2  seconds 
for  particle  sizes  of  5,  30  and  200  microns.  Even  with  the  leak,  the 
200-micron  particles  moved  in  a  slug-type  fashion,  as  shown  in  Figures  8e 
and  8f.  Figure  8d  shows  that  the  five-micron  powder  filled  the  complete 
bag  as  early  as  0.2  seconds.  No  voids  were  present  at  this  time  in  the 
computation,  in  contrast  to  that  for  the  standard  30-micron  powder. 

POWDER  UNIFORMITY 

A  computer  run  was  made  for  configuration  j,  with  two  opposing  jets 
at  the  top  and  for  the  standard  30-micron  powder.  Table  III  shows  the 
uniformity  that  can  be  achieved  with  this  size  powder  by  utilizing  such  a 
configuration. 

COMPARISON  TO  HIGH-SPEED  PHOTOGRAPHIC  RECORDS 

A  comparison  of  the  computed  profiles  to  those  observed  in  high-speed 

photographs  was  made  (See  Figure  9) .  The  dashed  curves  are  the  contours 
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of  a  solid  volume  fraction  of  2x10  .  This  is  ten  percent  of  the  incoming 

volume  fraction  of  the  solids.  At  a  time  up  to  0.040  seconds  these 
computed  contours  agreed  with  the  contours  obtained  from  the  movies  of 
the  experiment.  At  later  times  this  computed  powder  concentration  reached 
the  far  end  of  the  bag  earlier  than  shown  in  the  experimental  photographic 
records.  This  may  be  due  to  the  fact  that  in  the  experiment  the  bag  was 
partially  deflated  initially,  thus  offering  resistance  to  the  forward 
flow  of  powder.  After  powder  deflection  from  the  end  of  the  bag,  experi¬ 
mental  and  theoretical  results  again  are  in  agreement.  For  example,  the 


shape  of  the  computed  jet  closely  resembled  the  observed  jet  shape.  The 
theoretical  computations  and  the  experimental  results  both  showed  a 
narrowing  and  a  widening  of  the  jet. 


Figure  9  also  shows  the  computed  concentration  profiels  for  a  volume 
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fraction  of  the  solid  equal  to  18x10  at  0.040  and  0.070  seconds.  This 
is  90  percent  of  the  inlet  volume  fraction.  This  concentration  lags  con¬ 
siderably  behind  the  10  percent  of  inlet  volume  fraction.  Clearly,  the 
dispersion  of  the  powder  is  quite  large.  Hence  a  better  correlation  of 
theoretical  and  experimental  results  would  require  a  knowledge  of  the 
value  of  powder  concentration  at  the  visible  concentration  front. 

It  is  possible  that  the  computations  could  provide  an  estimate  of  the 
solid  volume  fraction  visible  in  high  speed  photographic  records.  In 
addition,  discrepancies  concerning  the  initial  shape  of  the  bag  and  the 
amount  and  location  of  leaking  exist  between  the  computer  modeling  and 
the  experimental  work.  As  these  discrepancies  are  resolved  better  agree¬ 
ment  between  the  theory  and  experiment  should  be  obtained. 
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NOTATION 


Cp  Drag  coefficient 

Diameter  of  solid  particles,  m 

2 

g  Gravitational  force  per  unit  mass,  m/sec 

G  Modulus  of  elasticity,  particle-particle  interaction 

coefficient,  N/m^ 

P  Pressure,  Pa 

Reg  Solids  Reynolds  number 

t  Time,  sec. 

U  Lateral  gas  velocity,  m/sec 

R 

Ug  Lateral  solid  velocity,  m/sec 

V  Axial  gas  velocity,  m/sec. 

S 

Vg  Axial  solid  velocity,  m/sec. 

x  Co-ordinate  in  lateral  direction,  m 

y  Co-ordinate  in  axial  direction,  m 

Greek  Letters 

8  and  8  Fluid-particle  friction  coefficient  in  the  x  and  y  directions 

x  y  ^ 

K  /m  -sec. 

g 

e  Gas  volume  fraction 

T  Solids  stress,  related  to  particle-particle 

pressure,  N/m*. 

y  Gas  viscosity.  Kg/ (m. sec) 

g 

3 

p  Gas  density.  Kg/m 

g 

3 

pg,Pp  Solid  and  particle  densities.  Kg/m 

<J>8  Sphericity 
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Table  II.  Solid  volume  fractions  for  configuration  "d"  after 

1.1  seconds  (top  table)  and  1.2  seconds  (bottom  table) 
(The  ables  show  the  values  at  the  left  half  of  the 
cav-'+  t  40cm  intervals  horizontally  and  6.1cm 
int  ..  ’’ertically) 
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Table  III.  Solid  volume  fractions  for  configuration  "J"  after  0.7 
seconds  (top  table)  and  1.0  second  (bottom  table)  (The 
tables  show  the  values  at  the  left  half  of  the  cavity  at 
40  cm  intervals  horizontally  and  6.1  cm  intervals 
vertically) . 
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400  cm 


(d)  Opposed,  Centered 


Figure  1A  Four  Beg  Dissemination  Configurations 
(Gas  velocity  =  40  ro/sec;  Particle 
diameter  =  30ym,  particle  density  =  1.65  gm/cc 
Inlet  solids  volume  fraction  =  0.002^. 
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400  cm 


88  cm 


Terminate  powder  at  0.82  sec 


(d)  Opposed,  Centered,  dfe 
No  Leak 


30.5  cm 


Leak 


terminate  powder  at  0.82  sec 


Opposed,  Centered, 
Uniform  Leak  across  Top 


30.5  cm 


Leak 


Terminate  powder  at  0.53  sec 


Centered,  Single, 
Uniform  Leak  across  Top 


Figure  IB 


Three  Bag  Dissemination  Configurations 
(Sinusoidal  gas  inlet  velocity  (m/sec)  = 

40  +  10.5  Sin(20irt),  where  t  is  in  seconds; 
particle  diameter  =  30ym;  particle 
density  *  1.65  gm/cc;  Inlet  solids  volume 
fraction  *  0.002) 
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Gas  Velocities  for  configuration  (e)  with 
a  leak  at  0.1  (top  graph),  0.2  (middle 
graph)  and  at  0.3  seconds  (bottom  graph). 
Scale:  1  cm  =  107  m/sec. 


DIST.  FROM  THE  TOP 

Fig.  6b.  Effect  of  Leak  on  Powder  Concentration  at  0.1  secon 
at  Various  Positions.  Legend:  l~no  leak,  configura¬ 
tion  d.  2  -  with  leak,  configuration  e. 


OIST.  FROM  THE  TOP 

Fig.  6e.  Effect  of  Leak  on  Powder  Concentration  at  0.1  seconds 

at  Various  Positions.  Legend:  1-no  leak,  configuration 
d.  2-with  leak,  configuration  e. 


DIST.  FROM  THE  TOPCCM) 

Fig.  7a.  Effect  cf  Leak  on  Powder  Concentration  at  0.2  seconds 
at  Various  Positons.  Legend:  1-no  leak,  configuration 
2-with  leak,  configuration  e. 
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Fig.  8e.  Powder  concentrations  for  200  ym  particles  at 
0.1  sec.  Configuration  i. 


Fig.  8f.  Powder  concentrations  for  200  ym  particles  at 
0.2  sec.  Configuration  i. 


Figure  8.  EFFECT  OF  PARTICLE  SIZE  ON  POWDER  VOLUME 
FRACTIONS . 


Figure  9.  Comparison  of  Computed  (configuration  k)  and  Experi¬ 
mental  (Run  BM-4)  Concentration  Profiles  at  Various 
Times  (Seconds)  after  Injection.  Solid  Line-Expeji- 
mental.  Dashed-Volume  Fraction  of  Solids  of  2x10 
Dot-Dash-Computed  Volume  Fraction  of  18x10"  . 
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